Charge ordering in half-doped Pr(Nd)o.5Cao.5Mn03 under magnetic field 
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Recent experiments in Pr(Nd)o.5Cao.5Mn03 thin films exhibited (multiple) reentrant charge or- 
dering (CO) transitions with change of temperature T under fixed magnetic field H, which are in 
contrast to the results for the corresponding bulk materials. To explain the experimental findings, a 
model including the double- exchange mechanism, intersite Coulomb interaction and electron-phonon 
coupling is proposed, for which the reentrant CO is naturally obtained due to the temperature de- 
pendence of the band-type coherent polaron hopping. Various results for the CO in the {H, T) 
plane are extensively discussed. The theory is considered to be valid for both thin films and bulk 
materials. 



PACS numbers: 75.30.Vn, 71.30.-Hh, 71.38.-k, 75.70.Pa 
I. INTRODUCTION 



The colossal mag- 

netoresistance manganites Ri_a;Aa;Mn03 (R: rare-earth 
element like La, Pr, Nd; A: divalent alkali like Ca, Sr) 
and R„_„a;An-„a:Mn„03„+i {n = 1,2) have attracted 
great interest for both bulk materials and thin films in 
recent years.Ej Due to the complex interplaypoLcharge, or- 
bital, spin, and lattice degrees of freedom, intriguing 
phenomena have been discovered. Among them, charge 
ordering (CO) with alternating Mn'^+ and Mn'^+Jons in 
real space is one of the most remarkable findings.^ It has 
been observed for optimal doping x = 0.5 in various (R, 
A) compositions, where the concentrations of Mn'^+ and 
Mn^+ ions are equal, and was also found in tie doping 
region 0.3 < a; < 0.75, depending on R and A.eI The CO 
may be melted by external magnetic fields, which are 
favorable to a ferromagnetic (FM) metallic state. 

The CO iu, paanganite thin films has been studied 
very recently.LTtll Due to strain effect, different proper- 
ties from those in bulk materials are expected. Prellier 
et al. studied the CO instabilit3£.jjnder magnetic field 
H in Pro.sCao.sMnOa thin filmsLlu and found a spec- 
tacular reduction of the melting field He for CO under 
tensile sttain. This is ascribed to an enhanced electron 
itineracy.EI On the other hand, the critical He, as a, func- 
tion of temperature T, is found to be non-monotonous. 
From the phase diagram in the plane {H, T)j3 which is 
replotted in Fig. ^ by the solid line, a more interesting 
phenomenon, i.e., reentrant CO transitions with change 
of T under fixed H , can be clearly seen. For example, at 
H = 6T, the CO appears at about T = 210 K but later 
mehs at T = 100 K. The profile for i? = 3.5 T is even 
more remarkable: the CO appears at about T — 220 
K, and later melts at T = 60 K, but finally reappears 
at r = 20 K, i.e., multiple (twice) reentrant CO tran- 
sitions are expected through the whole temperature re- 
gion. A similar phase diagram with possible single reen- 
trant behsftaor was also reported for Ndo.sCao.sMnOa 
thin filmsjlil and is reproduced in Fig. |l| by the dotted 
line. In contrast, no obvious reentrant jGO-was found for 
the two corresponding bulk materials. t3ii3 Naturally, a 



theoretical investigation of the experimental findings is 
required. 
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FIG. 1. The experimental CO phase diagrams in 
the {H, T) plane for Pro.sCao.sMnOs (PCMO) and 
Ndo.sCao.sMnOa (NCMO) thin films, extracted from Refs. 
[8] and [11], respectively. Note that the critical field has 
been taken as the average of the values in the field-increasing 
and decreasing scans. 

It should be pointed out that the reentrant CO 
under magnetic field was clearly observed in bulk 
Ri_:ECaa;Mn03 (R=Pr, Nd) forj-cipc|inc_.awa?/ /rom 0.5, 
e.g., 0.3 < a; < 0.45 for R=Pr.EjMElllJ Although this 
issue is not fully understood either, it was argued by 
Khomskii et al. that phase separation into an electron- 
rich region (due to extra Mn^+ ions) and a CQjjegion, 
as that for x — 0.5, is a relevant mcchanism£3'tHl On 
the other hand, even in the half-doped {x — 1/2) case, 
(multiple) reentrant CO was found before in bi-lavered 
manganite LaSr2Mn207 under zero magnetic field.Ef — 
A mechanism for this scenario was proposed by us.E^ 
The idea is to ascribe the CO to the competition of elec- 
tronic kinetic energy and intersite Coulomb interaction, 
and simultaneously take electron-phonon (e-p) coupling 
into account, i.e., consider polaron ordering. Due to the 
inherent temperature-dependent band-type polaron mo- 
tion, the kinetic energy may be small compared to the 
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Coulomb energy at some temperatures (favoring CO), 
but dominates at other temperatures (favoring the ho- 
mogenous state) so that a reentrant behavior may ap- 
pear. Here, the same idea will be generalized to study the 
CO under magnetic field. The experimental phase dia- 
grams are interpreted and, more generally, various results 
for CO under magnetic field for half-doped manganites 
are discussed. 

Before we start to work out the details it deserves to 
further clarify the validity of the above proposed the- 
ory. Actually, we notice that there are several alterna- 
tive models to study CO (and concomitant orbital and 
magnetic orderings) in half-doped manganites currently. 
A natural model is to take the intersite Coulomb inter- 
action as the driving force-f or M CQj This modelling was 
adopted by many author^"Einl3'E3 and also used for a 
discussion of other raajterials, e.g., the Verwey t ra i ng ition 



m mag netite Fe304.E3 Another kind of modelr fri em- 
phasizes the role of CE-type antifcrromagnetic (AFM) 
spin order concorpitant with CO, which was found in e 
Ndo.5Sro.5MnO3.El As derived by van den Brink et al. 
the on-site Coulomb interaction alone between electrons 
in different orbitals will lead to CO. And the experimen- 
tal (7r,7r,0)-type CO, i.e., the checkerboard CO forms in 
the a;y-plane and stacks in z-direction, can be explained 
within the same picture. However, this model is not well 
suited for Pr(Nd)o.5Cao.5Mn03(_svliere the spins are not 
yet ordered when CO appears.EJilj Recently, a purely 
Jahn- Teller (JT) phononic model, even without di 
electron correlations, was proposed by Yunoki et al. 
By treating the e-p coupling adiabatically, they have con- 
structed the phase diagram of charge, orbital and mag- 
netic orderings with change of JT coupling strength and 
AFM superexchange of the core spins. This model seems 
not to be generally accepted at the moment .E3 With- 
out relying on peculiar features of concrete materials 
(e.g., magnetic structure), the intersite Coulomb inter- 
action, which is always present in reality, is regarded 
as the origin of CO as done by most authors. And 
on the other hand, another aspect of e-p coupling, i.e., 
polaron formation with dynamic phonons is addressed. 
Experimentally it has been verified that the CCLahmild 
be viewed as the formation of a polaron latticeEjtj A 
strong evidence is the isotope effect, where the CO tran- 
sition tenmerature T^o was changed by replacing -'^^O 
with -'^^O.Ea Therefore, as the simplest model to study 
CO in manganites, the e-p coupling is indispensable, be- 
sides the intersite Coulomb interaction. This improves 
the one pmposed by Khomskii where only the latter is 
includedo Indeed, the model with both e-p and inter- 
site Coulomb interactions has been sucoeasfully used to 
explain the experimental isotope-effect J£j renormaliza- 
tion of sound velocity around T_B3 and reentrant CO in 
LaSr2Mn2 07 mentioned aboveHj Also, it is believed to 
have grasped much of the physics of the CO transitions 
in Pr(Nd)o.5Cao.5Mn03, addressed in this paper. 



II. MODEL HAMILTONIAN AND TECHNICAL 
TREATMENTS 

On the basis of the considerations in the introductory 
section, we write the Hamiltonian as 

H = -t{cos{9/2)) ^(cjcj + h.c.) + n^n^ 

i i 

The first term is the double-exchange (DE) model, where 
t is the nearest- neighbor (n.n.) hopping and 9 is the angle 
between two n.n. core spins which is assumed uniform or 
an average on Oij . V is the intersite Coulomb interaction 
and Tii = clci is the electron number operator. Irre- 
spective of the details for various types of e-p coupling 
(e.g., JT, breathing, etc.), we have used a simplexauaite 
Holstein-typc model with dispersionless phononsoEJ as 
shown by the third term to describe an essential effect 
of e-p coupling. Then bi, bj are local phonon operators 
at site i, g is the e-p coupling constant and cu is the 
Einstein phonon frequency. The band is half-filled, i.e., 
one electron per two sites. Note that the effect of ex- 
ternal magnetic field H has been implied in the factor 
(cos(0/2)), which becomes larger with increasing H due 
to enhanced FM order (i.e., 9^0). 

In the following, we consider the Hamiltonian (^ in 
the two-dimensional (2D) case since the checkerboard GO. 
was found only in the basal plane as mentioned above.E3 
Moreover, we mainly focus on the qualitative results in 
this work, thus no sophisticated techniques will be elab- 
orated. The quantitative aspects of the results obtained 
should not be trusted literally, but they are qualitatively 
correct. Throughout the paper we set fi ~ ks ~ ^ and t 
is taken as the energy unit. 

First, the usual Lang-Firsov (LF) transformation is 
acted upon Hamiltonian (0), i.e., H = UHW with 

U = exp[-{g/u;)J2Mh~bl)] . 

i 

Then, thermally averaging on the phonon states, one ob- 
tains the following effective electronic Hamiltonian: 



ifcff = -i^{c\cj + h.c.) + V^', 



(2) 



with the renormalized hopping i = t{cos{6 /2))P{T). 
Here 

P(T) = exp[-(^)2coth(;^)] 

is the polaron narrowing factor and a constant 
— {g^ /oj)^^nf is ignored in Hamiltonian (||). Thus the 
effect of e-p coupling is contained exclusively in P{T). 
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FIG. 2. P(T) as a function of T for g = a; = 0.1 (solid line) . 
The dashed line schematically shows the result with acoustic 
phonons (see text). 

Note that P{T) has a strong temperature dependence, 
which is plotted in Fig. |^. Although the narrowing ef- 
fect is overestimated by the LF method for uj/t <^ 1, 
i.e., the absolute value of P{T) is too small, the T- 
dependent property is qualitatively correct, i.e., P{T) 
increases monotonically with decreasing T because of 
smaller and smaller phonon occupation number. Due to 
the same reason, P{T) is seen to be nearly flat at very 
low temperatures, about T/lo < 0.2, where the phonons 
are hardly excited. T.f acQus^tic phonons are considered, 
as in some references Jij'oEflcj a rapid variation of P{T) 
with T should be still possible up to T 0, as schemati- 
cally shown by the dashed line in Fig. |[ We emphasize 
here that the T-dependent P{T) is the essence of our the- 
ory, which may render the CO to melt with decreasing 
T by quickly enhancing the hopping t. In principle, the 
thermal average (cos(0/2)) is also T-dependent, but it 
is assumed to be quite weak compared to P{T). Actu- 
ally in the CO states in Fig. |^, there is no indication for 
noticeable changes of the spin order under fixed, fields, 
especially in the reentrant temperature region.E3 Thus 
the possible tiny change of (cos(0/2)) with T, even in 
the correct direction (i.e., increasing with decreasing T), 
is argued impossible for a melting of the CO, and thus 
not important. In what follows, (cos(0/2)) will be simply 
considered as a function of H, i.e., (cos(0/2)) = F{H). 
Finally, a minor point is that, for the current model (|^), 
no phonon-mediated intersite interaction has been de- 
rived. For other types of e-p coupling, such an extra 
interaction rnajiJpe obtained, leading to a shift of the 
parameter y£jo 

For the effective Hamiltonian (H), we make use of the 
Hartree mean- field (MF) theory to study the CO. For a 
bipartite lattice with sublattices A and B, it is assumed 
_ / 1/2 + a;, iek 
^ 1/2-a;, i e B 
viation from the uniform electron distribution. The final 
self-consistent equation for the order parameter ob- 
tained as follows (trival solution .t = is ignored)! 



{n^) = 



where x characterizes the de- 



2V tanh[2V(iz)2 



ivxr/T] 



TT 



/o 



Viizy + (Vx) 



K{VT^) dz , (3) 



where K is the complete elliptic integral of the first kind. 



III. RESULTS AND DISCUSSIONS 

The above equation can be easily evaluated numeri- 
cally. For manganites, the bare hopping i is ~ 0.2-0.5 
eV, and the intersite inteKiction V/t is ^ 0.1-0.2 due to a 
large dielectric constant .l3 The optical phonon frequency 
u! ranges ^ 10-70 meV. For example, a soft mode at w ~ 7 
meV was reported by Zhae et al. for Lai_a;Caa:Mn03 
{x = 0.25, 0.4) thin films.N We point out that the high 
frequency modes play no role for the melting of CO here 
because, as discussed above, they only contribute to a 
nearly flat P{T) in the temperature scale for the reen- 
trance to occur. 

The solution for V = 0.1, oj = 0.05 is displayed in the 
phase diagram Fig. ^ in the plane of T and F{H). The 
results for different g values are presented. The reen- 
trant CO is clearly seen within a finite region of F{H) 
for g = 0.04, 0.042, and 0.043, while not for g = 0.035 for 
any < F{H) < 1. For a careful analysis of the reentrant 
behavior, the value g = 0.042 is taken as an example and 
enlarged in Fig. I We fix F{H) = 0.8 as indicated by the 
arrow and check how the CO state evolves with decreas- 
ing T. At the highest energy scale, the entropic term 
dominates the free energy and the homogeneous state is 
stable. For temperatures of the order of V, where the 
coherent polaron motion is nearly fmpen (i.e., t —^ 0), 
a charge order transition takes placeo On decreasing T 
the band-type hopping t becomes larger, and moreover, it 
increases fast enough to dominate V so that the CO state 
is destabilized. Thus the second transition, where the CO 
melts into the homogeneous state (hs) takes place. With 
continuing decrease of T, the third transition occurs to 
recover the CO since the ground state must be charge 
ordered, which will be discussed in detail later. To sum- 
marize, in the whole temperature region, the system goes 
through a series of transitions: hs CO ^ hs ^ CO 
with decreasing T. The corresponding x vs. T is shown 
in the inset of Fig. || by the solid line, where the order 
parameter takes on alternating zero and non-zero values. 
Interestingly, the double reentrance of the CO at fixed 
F{H) or external field obtained here is a remarkable ex- 
perimental observation shown in Fig. |^. In addition, in 
the inset of Fig. ^, the curves a; vs. T for two examples 
without reentrance, i.e., F{H) — 0.7, 0.9, are also shown. 
An analysis of them is helpful. For F{H) = 0.7, as shown 
by the dashed line, x decreases with decreasing T (due 
to increasing t) in the intermediate region but does not 
approach zero. This means that the increase of t is not 
sufficient to melt the CO. The case for F{H) — 0.9 is 
different. Now the absolute value of i becomes relatively 
large, even at high temperatures. Although it is reduced 
at higher T, the simultaneously enhanced thermal fluc- 
tuation will prohibit CO. Thus the CO can not be sta- 
bilized at a relatively high temperature as the cases for 
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F{H) = 0.7 and 0.8. It remains unstable with decreasing 
T, and appears only at very low temperatures to realize 
the charge ordered ground state, as seen by the dotted 
line. 
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FIG. 3. The phase diagram in the plane of T and F{H) 
for V — 0.1, uj — 0.05. The solid curves from left to right 
correspond to g = 0.035, 0.04, 0.042, and 0.043, respectively. 
The dashed line shows a correction to the phase boundary of 
g = 0.035 (see text). 
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FIG. 4. The enlarged phase diagram for g = 0.042 in Fig. |[ 
The two dashed lines 1 and 2 show two possibilities for a 
correction to the phase boundary (see text). In the inset the 
order parameters x vs. T are plotted for F{H) = 0.7, 0.8, 
and 0.9. 

As for the ground state, it is always charge ordered 
for any V/t > because of the perfect nesting tyoperty 
at half- filling for the model Hamiltonian (^.l3'E2I This 
conclusion can be drawn from Eq. (^), which is mxjven 
to always give a solution of nonzero x at T — 0e3 On 
the other hand, it must be noted that the perfect nest- 
ing property may be easily broken in real materials, for 
example, by inclusion of the next n.n. hopping which is 
usually present. In situations away from perfect nesting, 
a natural conclusion for the ground state is that the CO 



only appears when V exceeds a critical value, or equiv- 
alently it vanishes if i becomes larger^ Correspondingly, 
in the phase diagram of Fig. ^, a critical value of F{H) 
(or a critical field He) is expected at T = to separate 
the CO and homogeneous states. With this in mind, 
the phase boundary may be reasonably corrected in two 
ways, 1 and 2, as shown by the dashed lines. For line 1, 
only the possibility of one reentrance is obtained. And 
for line 2, a double reentrance is still possible. These 
two possible phase diagrams are qualitatively consistent 
with the experimental findings in Fig. ^ At the same 
time, a similar correction is valid for the phase bound- 
ary of g — 0.035 in Fig. |^, as shown by the dashed line. 
This corrected phase boundary is the same as usually ob- 
served in half-doped manaaiiitps in the [H, T) plane, e.g, 
the bulk Pro.5Cao.5MnO3.0jllil For concrete calculations 
including all above corrections, a model beyond perfect 
nesting is needed, which is left for future studies. 

As mentioned before, the reentrant CO on decreasing 
T depends on the increasing strength of t and its relative 
magnitude to V . Under which conditions it takes place 
is complicated, relying on the concrete parameters, i.e., 
g, w, F{H) which control t, and the value of V. Some 
insight can be gained by comparing the phase diagrams 
with different values of g, w, and V . A few consequences 
seem clear: for relatively large F(H), increasing g helps 
the first CO transition to occur at a high temperature, 
and consequently favors the appearance of the reentrant 
behavior; a similar effect can be obtained by increasing 
V; large uj is unfavorable to the reentrant CO. Unfortu- 
nately, nOjmiantitative criterion can be straightforwardly 
extracted.E^ 

So far, we have not addressed the observation that the 
reentrant CO is exhibited in Pr(Nd)o.5Cao.5Mn03 thin 
films, while, it is not obvious for the corresponding bulk 
matcrials.lis The difference is understandable from the 
above results. For thin films, due to strain effects, the 
lattice parameters are modified (e.g., the in-plane param- 
eter is elongated by the tensile strain). Correspondingly, 
the model parameters will change in thin films compared 
to bulk materials O We speculate that the reentrant be- 
havior in thin films is due to the softened optical phonon 
modes. 

Finally, we comment on the function F(H), which was 
used throughout the results. Its explicit form is not 
known although its qualitative behavior is transparent. 
Above, spin order was fixed when the reduced Hamilto- 
nian (|^) was introduced. Strictly speaking, the effect of 
H should be studied interactively with the CO under the 
same Hamiltonian by inclusion of the superexchange Jaf 
of the core spins and the Zeeman energy. The explicit 
form of F(H) can be roughly estimated in the following 
phenomenological way. We consider a Hamiltonian of the 
form -t{cos{9/2)) ^^^^-^ (clc, +h.c.) + Jaf • S, - 

g^J■BH J2i J with S: core spin and g: Lande factor. 
Moreover, a uniform canted state from the competition 
between FM order due to the DE mechanism and AFM 
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superexchange is assumed in order to simulate the es- 
sential effect of the external field (i.e., favoring the FM 
order). In this way, a liuear relation between (cos(0/2)) 
and H may be obtainedo 



IV. CONCLUSION 

Based on a double-exchange model with intersite 
Coulomb interaction and e-p coupling, charge ordering 
has been studied under magnetic field H in the context 
of polaron formation. A melting of the CO on decreas- 
ing T under fixed H can be naturally obtained due to 
the T-dependent band-type polaron hopping. More gen- 
erally, various phase diagrams in the plane of T and 
F{H) with (single or double) or without reentrant be- 
havior, depending on the concrete model parameters, 
have been extensively discussed. The experimental find- 
ings on the CO under magnetic field in the half-doped 
Pr(Nd)o.5Cao.5Mn03 can be explained for both thin films 
and bulk samples. Future work will have to address 
a number of important issues which are still not well 
comprehended: the orbital degree of freedom, a devia- 
tion from (tt, tt) for the wavevector of the CO in thin 
films J3 anjd-a possible magnetic phase separation in the 
CO state.El 
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